Superfluid and supersolid phases of lattice bosons with ring-exchange interaction 
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We examine the superfluid phase of a hard-core boson model with nearest-neighbor exchange 
J and four-particle ring-exchange K at half-filling on the square lattice. At zero temperature we 
find that the superfluid in the pure-J model is quickly destroyed by the inclusion of negative-if 
ring-exchange interactions, favoring a state with a (tt, it) ordering wavevector. Minimization of the 
mean-field energy suggests that a supersolid state with coexisting superfluidity, charge-density wave, 
and valence-bond-like order is formed. We also study the behavior of the finite-T Kosterlitz-Thouless 
phase transition in the superfluid phase, by forcing the Nelson-Kosterlitz universal jump condition 
on the finite-T spin wave superfluid density. Away from the pure J point, Tkt decreases rapidly 
for negative K, while for positive K, Tkt reaches a maximum at some K ^ in agreement with 
recent quantum Monte Carlo simulations. 



I. INTRODUCTION 

The road to discovering new quantum phenomena 
such as exotic phases or unconventional quantum phase 
transitions is paved by the understanding of quantum 
Hamiltonians with competing microscopic interactions. 
Recent work has shown that competing kinetic energy 
terms, particularly four-particle "ring-exchange" , is a fer- 
tile venue for the study of valence-bond solids and uncon- 
ventional quantum phase transitions*^^ 

Long relegated to toy models of easy-plane spin sys- 
tems or simple boson theories, the understanding of U(l) 
Hamiltonians has resurfaced as critical for the engineer- 
ing of exotic quantum phases of ultracold atoms in op- 
tical lattices. Buchler el presented the design of a 
ring-exchange interaction for bosonic cold atoms in two- 
dimensional (2D) square optical lattices. There, a con- 
ventional nearest-neighbor hopping term competes with 
a four-particle ring-exchange, which arises out of the hop- 
ping of boson pairs on the corners of square lattice pla- 
quettes to an intermediate molecular state and back to 
the (oppositeVplaquette corners. The Hamiltonian pro- 
posed in Ref. |5j can be written, 



H 



-J 



-K 



(ijki) 



where the first term is the usual boson hopping, and the 
second is the ring-exchange interaction. Here, i and j are 
neighboring sites, lying opposite to sites I and k (respec- 
tively), which together form the basic square plaquette 
of the 2D lattice. This model has drawn significant at- 
tention recently from proposals that it may harbor an 
exotic quantum liquid phase that possesses d-wave cor- 
relations - a so-called d-wave Bose liquid (DBL) phase 6 
- including extensive numerical investigation using exact 
diagonalization and the density matrix renormalization 
group Z 

The main difficulty with the Hamiltonian Eq. dXJ) that 
prevents an exact solution in 2D for large systems (or the 
thermodynamic limit) is the presence of the prohibitive 
"sign-problem" in quantum Monte Carlo (QMC) simula- 



tions for the parameter regime K < 0. In fact, the same 
model in the parameter regime without the sign prob- 
lem, K > 0, has been solved by QMC previously both at 
half-filling 2 and in the presence of a symmetry-breaking 
chemical potential^ In this model, it was demonstrated 2 - 
that the superfluid phase that dominates for large J is de- 
stroyed by increasing the magnitude of the ring-exchange 
interaction, which realizes first a (tt, 0) valence-bond solid 
(VBS) for 8 < K/J < 14, and a (tt,tt) charge-density- 
wave (CDW) for K > 14. The intermediate superfluid to 
VBS quantum phase transition was studied intensely 10 
as one of the first candidates for a deconfined quantum 
critical point. 9 However, other interesting behavior oc- 
curs in the model, in particular in the experimentally 
relevant regime of finite-temperatures, where relatively 
little attention have been paid^ 

In order to elucidate the mechanism by which the su- 
perfluid phase, which occurs for dominant J in Eq. |T]), 
is destroyed by the four-site K interaction, we explore in 
this paper the behavior of the model using linear spin- 
wave (SW) theory. First, we calculate the dispersion and 
superfluid density p s at T = 0, where its behavior as a 
function of K indicates the realization of a phase with 
(n, 7r) symmetry for K/J < —2. Exploring the order- 
ing nature of this phase in mean-field theory suggests 
a (tt, tt) modulated in-plane order parameter, the mean- 
field indicator of a supersolid state, where superfluidity 
is expected to coexist with both CDW and also valence- 
bond-like order. In the superfluid phase, we also calcu- 
late the superfluid density at finite-temperature, and es- 
timate the position of the Kosterlitz-Thouless (KT) tran- 
sition using the universal jump condition. Interestingly, 
for K > the linear spin-wave theory reproduces the 
QMC result^ - that the maximum temperature of the KT 
transition does not occur when K = 0, rather it happens 
at some intermediate value of K/J. 



II. LINEAR SPIN WAVE THEORY AT T = 

In order to study the destruction of the superfluid 
phase by the competing ring-exchange interaction param- 



2 



eterized by K, we begin by examining the behavior of the 
system at T = in a linear spin- wave analysis. To begin, 
it is helpful to write our hard-core boson Hamiltonian as 
the analogous two-dimensional spin- 1/2 XY model using 
the standard mapping: b\ -> S+ and h -> Sf. This 
transforms the Hamiltonian Eq. (JTJ) into 

H=-J^2(SfSr + SrS+) 

-if £ (SfS]S+Sf + Sr S +S-S+), (2) 

which can be written with and 5* y operators as 



H 



2(7 / 



SfSf) 



(3) 
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As mentioned above, the groundstate of the pure-J 
Hamiltonian is a bosonic superfluid, or an in-plane fer- 
romagnet in the spin language, with an order parameter 
(S x ) at zero temperature. We therefore perform our 
spin-wave expansion around this ordered state, treating 
the K term of the Hamiltonian as a perturbation. As first 
demonstrated by Gomez-Santos and Joannopoulos^i the 
proper Holstein-Primakoff representation in the case of 
the XY model is 



Sf 

sf 



1 



a]a,i 



^(at-a,-), 



(4) 



giving us a leading order approximation to the Hamilto- 
nian in terms of bosonic spin- wave operators dj and a\ . 
Writing the Hamiltonian in Fourier space gives the linear 
spin-wave Hamiltonian 



are defined in terms of 



Y k 



2"4 7k ' 



- cos k x cos k v 
4 



where 



1 1 

2 
1 

2 7k 

7 7k - -. cosk x cos ky 
4 4 y 



7k = cos k x + cos k y 



(9) 
(10) 

(11) 
(12) 

(13) 



This form for the SW Hamiltonian reduces to that ob- 
tained by Bernardet et. alJ£ in the limit of the simple XY 
model (K = 0) and the absence of an external magnetic 
field. Following Ref. [IH, we diagonalize Eq. §5§ using a 
Bogoliubov transformation, 



(14) 



where ak and are destruction and creation operators 
for quasiparticles of momentum k. The bosonic commu- 
tation relations are satisfied if 



Uk = cosh(y> k ), v k = sinh(</? k ) 



(15) 



where ip k is determined by the requirement that it diag- 
onalizes our Hamiltonian. This yields 



k 2 



2 1 

k 2 



- 1 



and thus the diagonalized spin-wave Hamiltonian 12 
H 



H 



MF 



A, 



(16) 
(17) 

(18) 



Ev A k- B k(4"k + aL k a-k). 



H = H M f + E [ Ak ( a ic a k + aL k a-k) 



B k {ala^y. + a k a_ k ) . (5) 



It is the diagonal SW Hamiltonian, in this generalized 
form, that we now proceed to use to analyze the behavior 
of our superfluid phase at zero temperature. 



Here, the mean-field energy is 



A. Dispersion 



H 



MF 



-JN 



KN 



(6) 



We begin by studying the dispersion, which is obtained 
immediately from Eq. (fTS")) 



while the coefficients A k and £>k 

A k = JU k + KV k 
B k = JWk + KX^ 



2JAI 



(19) 



(7) This can be examined as a function of K/ J for the Hamil- 

(8) tonian Eq. For K — and k — > 0, one reproduces the 
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FIG. 1: (Color online) The dispersion wt, as a function of 
k = (k x ,k y ), for K/J = 2 (top) and K/J = -2 (bottom). 



expected linear dispersion^ a feature which survives for 
moderate K > (at least to first-order in the SW the- 
ory) - see Fig. [TJ No soft modes develop in the K > 
dispersion until very large values of K/J w 10 3 , which 
is well above the critical value of K/J = 7.91 where a 
phase transition to a (tt,0) VBS is known from QMC. 2 
In contrast is the behavior of the dispersion in the K < 
regime, which is intractable to QMC methods due to the 
negative sign problem. In this case, the dispersion reveals 
the development of soft modes at k = (tt,tt). The value 
of tends towards as K/J approaches —2, as il- 

lustrated in Fig. [1] This indicates that an ordered phase 
with (7r, 7r) symmetry is realized for sufficiently large \K\. 
The nature of this ordering is discussed in Section IIIII 



B. The zero-temperature superfluid density 

The most important quantity characterizing the super- 
fluid phase is the superfluid density (or spin stiffness), 
defined as the second derivative of the free energy with 
respect to a uniform twist (f> of the in-plane spin com- 



ponents across the system. At zero-temperature, one re- 
places the free energy by the groundstate energy, hence 
we use the definition 



d 2 E(<p) 

d 2 4> ' 



(20) 



which is valid for (f> — > 0. This is equivalent to finding 
the incremental energy per spin resulting from the twist, 



AE 
~N~ 



(H(cf>)) (H(Q)) 



N 



N 



1 

= ~ 2 Ps 



(21) 



Thus, we are faced with the task of finding H{4>). To do 
this, we consider the spins in a site-dependent rotated ref- 
erence frame, defined via the standard rotation operator 
about the z axis. 13 The relevant transformations become 



SJ cos 4>j 
Sj sin <f>j ■ 



■ Sj sin <fij , 
Sj cos 4>j , 



(22) 



which can also be written in terms of the spin raising and 
lowering operators, 



(23) 



We assume a uniform twist across the spins in the system, 



so that the twist 



= cf>, where i and j are nearest 



neighbors in the x or y direction. We find that, using 
the labeling of Eq. ([T]) , = (p^, fa = 0/ , and thus the 
<f> dependence of the ring-exchange term cancels. This 
gives for our twisted Hamiltonian 

H(<j>) = [(SfS* + S?SV) cos0 



+(SfSy -SfSf) sin <£\, 

-K 2_j (^i Sj S^S [ + S t Sj~S k Si~). 
(ijkl) 



(24) 



The linear components of the term proportional to sin cj) 
cancel, so in linear spin wave theory the twist effectively 
scales the nearest neighbor exchange value J — > J cos <j>. 
Thus, H(4>) is obtained directly from the Hamiltonian 
Eq. (0, with the redefinition of 



Hmf(<I>) = — JNcos( 

^4k(0) = JUfr COS (j) 

Bk(4>) = JW k cos 4 



KN 

~ "IT' 

KV k , 
-KX k , 



(25) 

(26) 
(27) 



but with the coefficients £/k, Vk, Wk, and X k unchanged 
from Eqs. (|9"1 fT2"|) . Expanding cos as 1 — 2 /2, and noting 
that at zero temperature, the thermal expectation value 
is 



n k 



(a k a k ) 



(al k a_ k ) = 0, 



(28) 
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FIG. 2: (Color online) The superfluid density as a function of 
K/J. The dashed line is the mean-field result, p s = 0.25. The 
linear SW result for K/J = is p a — 0.2709, which can be 
compared to the best known exact results p a = 0.2696(2). 14 



1 > — 

w 


k » 


> 
















^1 ' 




^1 


( i 
< 


I ^ 

<£1 


< — 



FIG. 3: (Color online) The (it, n) supersolid phase as dictated 
by mean-field theory, for K/J < —2. Solid (blue) arrows 
show the angles of deviation of the classical spin vectors from 
the x direction (the out-of plane deviation is zero). Dashed 
(green) arrows are the corresponding bosonic charge currents 
J, Eq IP), 



we find that 

* = i + 2]vE[ JC/ *-(^-^r 1/2 ( 29 ) 

k 

• (J 2 (Ul-W£) + JK(U k V k -W k X k )) . 

Note that we have divided by a factor of two to account 
for the fact that Eq. (fT8| is defined with respect to a 
twist in a single lattice bond, while Eq. |22|) considers a 
twist in each bond in the lattice, a total of 2N bonds. 

The T — superfluid density calculated in this lin- 
ear SW theory, Eq. (f29|) . is plotted in Fig. [2] for a range 
of K/J. Note that, as in the above section, we have 
set J = 1/2 to correspond to the usual definition of the 
XY spin model (when K — 0). The superfluid density 
curve has its maximum at K/J — 0, with a value of 
p s = 0.2709, which can be compared to the best nu- 
merical estimate for the groundstate spin stiffness in the 
XY model using finite-size scaling quantum Monte Carlo 
techniques, p s — 0.2696(2)^ Away from this maximum 
at K — 0, the superfluid density declines monotonically 
with increasing \K\. On the positive K side, p s decreases 
relatively gradually, only becoming zero for an extremely 
large value of K « 10 3 . This is consistent with the above 
results for the dispersion, which indicate that no soft 
modes develop for moderate values of positive K. Again, 
exact QMC results have revealed a T = quantum phase 
transition to a valence-bond solid phase at K/J « 7.91,- 
reiterating the limitations of the linear SW theory calcu- 
lation in this regime. 

On the negative-if side, where no QMC results are 
available, the value of p s drops rapidly as it approaches 
-2. In SW theory, there is a divergent negative contribu- 
tion to p s at this value of K/J = —2. This is of course 
an artifact of linear SW theory. As explained in the next 
section, one actually expects p s to remain finite on the 
other side of the transition, where an ordered phase with 
(71", 7r) symmetry is realized. 



III. SUPERSOLID PHASE FOR K < 

Examining our linear SW theory results of the previous 
section, for K/J < — 2 we expect the development of a 
phase with (tt,tt) as the ordering wavevector. We can 
further elucidate the nature of the ordering in this phase 
by minimizing the mean- field energy of the ring-exchange 
Hamiltonian, Eq. fT]). under the constraint of order with 
(■7T, 7r) symmetry. By representing the quantum spins as 
classical spin vectors, 

Sf = isin^cos^, (30) 
Sf = - sin 0i sin (f> i: (31) 
St = icos^, (32) 

we rewrite Eq. ([3]), which simplifies to the following 
Hamiltonian: 

H = — — ^2 sin^i sin6>j cos(0i — <fij) (33) 

— — \ s'mOi sin8j sin 9k sin 8i cos(</>,; — 4>j + 4>k — <Pi)- 

(ijki) 

Note that the angle of the classical vector from the x axis, 
4>i, should not be confused with the phase twist variable 
defined in the previous section. Restricting ourselves to 
half-filling, we look for solutions to the angles 0, and 4>i 
in the form 

9i = 9 + 9 ie ™ r <, (34) 
<t>i = 0o + 0ie l7r - r ', (35) 
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(where ir = (71", 7r)) by minimizing the classical energy 
Eq (J33J). The solution is 



% = e 1 = 00 = 

h = iarocos^V 



(36) 
(37) 



This solution corresponds to a staggered canting of spins 
away from the x axis, as illustrated in Fig. [3] In the boson 
language, we can map this to a charge current using the 
definition 



Jsin( 



This state would correspond to a staggered charge cur- 
rent as shown in Fig. [3] 

In MFT, the above state does not contain any charge 
density modulation, since the out-of-plane angle 9% is 
found to be zero. However, one expects (e.g. based 
on Landau-Ginzburg type symmetry arguments) such 
charge density modulation to appear once fluctuations 
are included. Such a state is therefore a bosonic su- 
persolid phase, with coexistent superfluidity and (71", w) 
CDW order. Interestingly, from the modulation of the 
boson current operator J illustrated in Fig. [31 it is clear 
that the supersolid phase in this case has the curious 
additional feature of valence-bond-like order coexisting 
with CDW density modulation, a situation which does 
not typically occur in lattice supersolids arising from such 
simple Hamiltonians. It is also worth noting that the 
mechanism for formation of this supersolid is distinct 
from that of other hard-core boson supersolid phases 
at commensurate filling discussed previously ; 1 5 1 1 6 1 1 7 i 1 8 
which are induced by the geometrical frustration of the 
lattice. 



IV. FINITE-TEMPERATURE 
KOSTERLITZ-THOULESS TRANSITION 



We turn now to a discussion of the supcrfluid phase in 
the J-K model at finite temperatures. The goal of this 
section will be to map out the finite-temperature phase 
boundaries of the superfluid phase in our linear SW the- 
ory, to compare (at least in part) to the results obtained 
in other studies. In the QMC work of Ref. H3, the phase 
boundary reported has a non-trivial shape, reaching a 
maximum for a positive value of K ^ 0. We are in- 
terested in whether this feature can be captured by our 
SW theory, therefore elucidating the physical mechanism 
by which this maximum at K ^ occurs. To do this, 
we first develop an expression for the superfluid density 
for T > (Section IIV A[) . then use it to estimate the 
SW Kosterlitz-Thouless transition through the universal 
jump condition (Section IIV Bp . 



A. The finite-temperature superfluid density 

At finite temperatures, we may calculate the superfluid 
density by recalling its original definition as the response 
of the free energy with respect to a twist. Calculating 
the partition function from Eq. (|18p we obtain the free 
energy 



f = h mf ^) + (A(^) 2 - s k (0) 2 - M4>)] 

k 

+ T]T hi (l-e-^/ T ), (39) 



/). (38) with the dispersion redefined as 



w k (0) = 2^ k (0) 2 - S k (0) 2 



(40) 



As noted previously, the twist-dependent terms in -ff (</>) 
are proportional to cos 4>, so we can write for (f> — * 



2N 

F = F(<j) = 0) + — p s <f 



(41) 



which is properly normalized with 2N being the number 
of lattice bonds. Therefore, we obtain for the superfluid 
density 



Pi(T) = lim — — — ^ 



dej) 2 



/T 1 duj^) 



(42) 



(43) 



where the important limits are evaluated as follows: 



lim 


9A k (0) 


--4 




4>^o 


d<j> 2 






lim 




7 ^k 






d(f> 2 


~ 2 




lim 


duj k ((j)) 


4 / 

= lim — A k 


dA k (cf>) 


<t>-*0 


dej) 2 


= -2 J— C/ k 4 
o>k 


dcj> 2 

S k 
-2 J— 1 



(44) 
(45) 



(46) 



Thus, this calculation yields the following expression for 
the superfluid density at finite-T: 



k 



2 2N 



(47) 



This SW expression for p s as a function of temperature 
is plotted for several parameters K/Jva Fig.|U 
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FIG. 4: (Color online) The superfluid density, calculated from 
Eq. (|47[1 . as a function of the temperature. Lines are labeled 
by the parameter value Kj J. The universal jump (Eq. (|48l> ) is 
illustrated as a dashed line. The exact value for the Kosterlitz- 
Thouless transition for K/J = is T KT = 0.6860(2), from 
unbiased QMC calculations.— 



B. Estimation of the Kosterlitz-Thouless transition 

The superfluid density of Eq. (|4"T)) . plotted as a func- 
tion of T, decays slowly (see Fig. [4]), crossing zero for rela- 
tively large temperatures. However, as is well known, the 
XY model realizes a Kosterlitz-Thouless^ (KT) transi- 
tion in any two-dimensional lattice, which is manifest as a 
discontinuity in p s precisely at Tkt- This, of course, can 
not be captured by a simple spin wave theory. We can, 
however, take advantage of the fact that the expected 
discontinuity in p s in 2D models such as this obeys the 
so-called universal jump condition, 

Tkt = \ps{Tkt), (48) 

first found by Nelson and Kosterlitzi^i One can thus 
hypothesize a reasonable estimate of Tkt by solving 
T = irp s (T)/2, using p s (T) from our spin- wave theory. 
An important test of this idea can be performed at K = 
and J = 1/2, the parameter values for the XY model. 
This is illustrated in Fig. [3J where the crossing point of 
the dashed line and the curve for K/J = is our SW 
theory solution to Eq. (|4"5|) . The value of Tkt obtained 
by this procedure is 0.6507, remarkably close to the exact 
result of T KT = 0.6860(2), obtained from unbiased QMC 
calculations 

We thus generalize this procedure of calculating Tkt 
in our linear SW theory to nonzero values of K/J. Re- 
sults are plotted in Fig. [5l As is evident by studying 
Fig. d] closely, the SW theory reproduces the remarkable 
trend that Tkt increases for small to moderate values 
of K > 0. In Fig. [4] the maximum in Tkt occurs at 
K/J w 6.5, before beginning to drop slowly. The phase 
boundary for the KT transition drops to zero for the large 
value of K/J « 10 3 , well outside of the expected range 




-2 2 4 6 8 10 
K/J 

FIG. 5: (Color online) The KT transition phase boundary 
as calculated in linear SW theory. A maximum of Tkt /J~ 
0.751 occurs at approximately K/J = 6.5. 

of validity for the SW calculation. Clearly, a simple spin- 
wave theory cannot capture the physics of the superfluid- 
VBS quantum phase transition that is observed in QMC 
simulations. 

On the negative- if side, the phase boundary drops 
rapidly to zero as K/J approaches —2, confirming the 
expectation from T ~ SW theory that a transition to 
a (ir, 7r) supersolid phase, discussed in the previous sec- 
tions, happens at this point. 

V. DISCUSSION 

We have studied through linear spin wave theory the 
various mechanisms which destroy the superfluid phase 
in the bosonic ring-exchange model Eq. (Q]) at half-filling, 
motivated as a candidate to realize exotic phases in cold 
atomic gases in optical lattices^ As is known from previ- 
ous QMC results,— sufficiently large ring-exchange K > 
promotes: first, a quantum phase transition to a valence 
bond solid state; second, a quantum phase transition to 
a (7r,7r) charge-density wave (CDW) state, at large K. 
In this paper, we have shown that a moderate value of 
K < is also sufficient to destroy the superfluid phase, 
promoting a (7r, 7t) ordered phase that is identified as 
a supersolid, or a phase with coexisting superfluidity, 
CDW, and valence-bond-like order. 

It is interesting to combine the exact results from QMC 
and our current SW theory results to speculate on the 
T = phase diagram for all parameters K/J at half- 
filling. The supersolid phase that we find at K/J < — 2 
must be related to the CDW phase found by exact QMC 
simulations for K/J > 14, since in the limit K/J — > oo, 
the Hamiltonian Eq. (JTJ) for K < maps exactly onto 
the Hamiltonian with K > 0, via a rotation of spins on 
one of the sub lattices by 7r/2 around the z axis. For 
a finite value of K/J < 0, two possibilities exist for the 
evolution of the supersolid state to the CDW; the state 
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can remain supersolid for all finite K/ J, or there could be 
a quantum phase transition between the supersolid and 
the (tt, tt) CDW solid at some critical Kj J. In either 
case, we have identified in this paper that a supersolid 
phase is possible at half-filling in the relatively simple 
Hamiltonian Eq. |T]) that has already been identified as 
a candidate for construction in ultracold bosons in optical 
lattices^ 

Further, we have studied the finite-temperature phase 
boundary of the superfluid phase in linear SW theory, 
by calculating the superfluid density and estimating the 
Kosterlitz-Thouless transition temperature by forcing it 
to obey the universal jump condition. We find that this 
procedure yields good numerical agreement with exact 
QMC results for K = 0. For K < 0, the phase bound- 
ary monotonically decreases to T — at K/J = —2, 
indicating promotion of the (71", tt) supersolid phase dis- 
cussed above. Increasing K > from the XY-point, 
the phase boundary initially increases, reaching a max- 



imum for K 7^ before monotonically decreasing. This 
is in qualitative agreement with the trend identified with 
QMC simulations, 10 and allows us to attribute at least 
the initial increase in Tkt for small K/J to the physics 
of non-interacting spin waves. It is remarkable that such 
nontrivial information about this model, namely the non- 
monotonic behavior of Tkt as a function of K/J, and the 
existence of a novel supersolid phase, can be contained 
in such a simple analytical theory. 
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